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Resummed ^(ajL n ),ff(a s aL n '),ff(a 2 L n ") corrections for n = 0, 1,2, n' = 0, 1,2, ra" = 1,2, 
in the presence of parton showers, on an event-by-event basis, without double counting and with 
exact phase space are required [jj], ® for precision LHC physics(l% or better total theoretical pre- 
cision [J3j] ) - We present here the first step in realizing our new MC event generator approach to such 
physics with amplitude-based QED^QCD resummation [Q] in the HERWIG6.5 [||] environment 
by our new parton shower MC for QCD, HERWIRI 1.0(31) [g]. HERWIRI 1.0(31) already shows 
improvement in comparison with the FNAL rapidity and soft pj data on single Z production as we 
quantify below. While the explicit IR cut-offs in the HERWIG6.5 environment will not be removed 
here, HERWIRI only involves integrable distributions so that these cut-offs could be removed. 

We first review our approach to resummation, which can be shown [J to be equivalent to 
those in Refs. [Q, 0], before we turn to a summary of the attendant new IR-improved DGLAP- 
CS [10, 11 1 theory [0] and a description of the implementation of the new IR-improved kernels in 



the framework of HERWIG6.5 [@] to arrive at HERWIRI 1.0(31). We illustrate the effects of the 
IR-improvement and compare with recent data from FNAL 1 . 

In Refs. [||, ^|] we have derived the following expression for the hard cross sections in the SM 
SU 2L x Ui x SU$ EW-QCD theory 

. - _ sum ir (qced) v 1 f d3 P 2 d3 i 2 ft d3kji ft ^ 3fc/ ^ 

n,m=0 n - m - J Pi 12 ji=l K h y 2 =l K J2 



(IK)* 



where the new YFS-style [ 13] residuals j3„ m (fci , . . . ,k n ;k[ , . . . ,k' m ) have n hard gluons and m hard 



photons and we illustrate the generic 2f final state with momenta p2 , qi specified for definiteness. 
The infrared functions SUMir(QCED), D QCE d are defined in Refs. [@, Eq. (|) is exact to all 
orders in a and in a s . 

The result Eq. ([]]) allows us to improve [f7|] in the IR regime the DGLAP-CS [10, 11] kernels 
as follows, using a standard notation: 



P% p (z)=C F F YFS ( Yg )el 5 * 



i—z 



P ^ iz)= C F F YFS ( 7q )e^ 1 + {1 Z) ' 



Gq 

P^( z )^2C G F YFS ( 7G )e^{ l —^ z ro + _£_(! _ ; 



+ ^(z 1+rc (l -z)+z(l -z) l+ K)-f G {y G )8{\ -z)}, 

P e q o (z) = FY FS (7G)e^ l -{z 2 (l -z)K + (1 -z)V G }, (2) 

where the superscript "exp" indicates that the kernel has been resummed as predicted by Eq. (|I|) 
when it is restricted to QCD alone and where we refer the reader to Refs. [^] for the detailed 



'FromRef. [|l|] the current state-of-the-art theoretical precision tag on single Z production at the LHC at 14 TeV 
is (4.91 ±0.38)% = (2.45±0.73)%(QCD + EW) ®4.11%(PDF) © \.\0±0M%(QCDScale) and for single W + (W") 
5.05±0.58%(5.24±0%). 
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definitions of the respective resummation functions Fyfs, 7a, 8a,/g,A = 1, G 1 . See Refs. [Q, [7|] for 
discussion of illustrative results and implications of the new kernels and Eq. ([I]) that are beyond the 
scope we have here. 

We have implemented the new IR-improved kernels in the HERWIG6.5 environment to pro- 
duce a new MC, HERWIRI1.0, which stands for "high energy radiation with IR improvement" 3 . 
We modify the kernels in the HERWIG6.5 module HWBRAN and in the attendant related mod- 
ules 4 with the following substitutions: DGLAP-CS P AB IR-I DGLAP-CS while leaving the 
hard processes alone for the moment. We have in progress [p^the inclusion in our framework 



of YFS synthesized electroweak modules from Refs. qi7|]for HERWIG6.5, HERWIG++ [|18p and 
MC@NLO [Il9|] hard processes 5 , as the CTEQ [|1 and MRST(MSTW) [§3|] best (after 2007) par- 



ton densities do not include the precision electroweak higher order corrections that are needed in 
a 1% precison tag budget for processes such as single heavy gauge boson production in the LHC 
environment 

The details of the implementation are given in Refs. [Q, ||] and we do not reproduce them here 
due to a lack of space. We have done many comparisons of the properties of the parton showers 
from HERWIG6.510 and HERWIRI1.031. In general, the IR-improved showers tend to be softer 
in the energy fraction variable z = E/Eseam where E {E^eam) is the cms parton(beam) energy for 
hadron-hadron scattering respectively. See Refs. [|], ||] for the complete discussion of such 
comparisons. We show in Fig. IT] [0] comparison analyses with the data from FNAL on the Z 



rapidity and pj spectra as reported in Refs. [ ]24] , |25| ]. We see that HERWIRI 1.0(31) and 
HERWIG6.5 both give a reasonable overall representation of the CDF rapidity data but that 
HERWIRI 1.031 is somewhat closer to the data for small values of F(The # 2 /d.o.f is 1.77(1.54) for 
HERWIG6.5 (HERWIRI1.0(31)).). Including the NLO contributions to the hard process via 
MC@NLO/HERWIG6.510 and MC@NLO/HERWIRI 1.031 [III 6 improves the agreement for 



both HERWIG6.510 and for HERWIRI1.031 (the # 2 /d.o.f are changed to 1.40 and 1.42 
respectively). That they are both consistent with one another and within 10% of the data in the 
low Y region is fully consistent with expectations and is an important cross-check on our work. A 
more precise discussion at the NNLO level with DGLAP-CS IR-improvement and a more 



complete discussion of the errors will appear pq]. We also see that HERWIRI 1.031 gives a better 
fit to the DO p T data compared to HERWIG6.510 for low p T , (for p T < 12.5GeV, the # 2 /d.o.f. are 
~ 2.5 and 3.3 respectively if we add the statistical and systematic errors), showing that the 
IR-improvement makes a better representation of QCD in the soft regime for a given fixed order in 



perturbation theory. Adding the ff(a !i ) correction from MC@NLO | |19| ] improves the # 2 /d.o.f for 
the HERWIRI1.031 in both the soft and hard regimes and it improves the HERWIG6.5 10 £ 2 /d.o.f 
for pj near 3.75 GeV where the distribution peaks. For pj < 7.5GeV the ^ 2 /d.o.f for the 
MC@NLO/HERWIRI1.031 is 1.5 whereas that for MC@NLO/HERWIG6.510 is worse. We 
await further tests of the new approach, both at FNAL and at LHC. - One of us (B.F.L. W) 
acknowledges helpful discussions with Prof. Bryan Webber and Prof. M. Seymour and with Prof. 
S. Frixione. B.F.L. Ward also thanks Prof. L. Alvarez-Gaume and Prof. W. Hollik for the support 



2 The improvement in Eq. (g) should be distinguished from the resummation in parton density evolution for the 
"z — > 0" Regge regime - see for example Ref. [ |l4[ |l5| . This latter improvement must also be taken into account for 
precision LHC predictions. 

3 We thank M. Seymour and B. Webber for helpful discussion on this point. 

4 We thank M. Seymour and B. Webber for helpful discussion. 

5 Similar results for PYTHIA | ^o| | and for the new kernel evolution in Ref. [ pT| ] are under study. 
6 We thank S. Frixione for helpful discussions with this implementation. 
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Figure 1: Comparison with FNAL data: (a), CDF rapidity data on (Z/y*) production to e + e pairs, the 
circular dots are the data, the green(blue) lines are HERWIG6.510(HERWIRI1.031); (b), DO p T spectrum 
data on (Z/y*) production to e + e~ pairs, the circular dots are the data, the blue triangles are HERWIRI1.031, 
the green triangles are HERWIG6.510- in both (a) and (b) the blue squares are MC@NLO/HERWIRI 1.031, 
and the green squares are MC@NLO/HERWIG6.510. These are untuned theoretical results. 



and kind hospitality of the CERN TH Division and of the Werner-Heisenberg Institut, MPI, 
Munich, respectively, while this work was in progress. S. Yost acknowledges the hospitality and 
support of Princeton University and the CERN TH Division. 
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